Transitions of two baryons to the H dibaryon in nuclei 
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We calculate the suppression in the rate at which two baryons in a nucleus (viz., nucleons or A's) 
convert to an H dibaryon, using an Isgur-Karl wavefunction for quarks in the baryons and H, and a 
Bethe-Goldstone wavefunction for the baryons in the nucleus. If rjj < 1/3 rjv, we find t Aaa _ a i h > ta 
and the observation of A decays from double-A hypernuclei does not exclude the existence of the 
H. If mu < 2m p , nuclei are unstable but have very long lifetimes. For reasonable values of rn 
and the nuclear wavefunction, the lifetime can be long enough to evade anticipated SuperK limits 
ta nn —,a' h >fewl0 29 yr, or short enough to be observed. An analysis of SuperK data to look for 
this possibility should be undertaken. 



I. INTRODUCTION 

The H dibaryon corresponds to the most symmetric 
color-spin representation of six quarks (uuddss). It is 
a flavor singlet state with charge 0, strangeness -2 and 
spin-isospin-parity I(J P ) = 0(0 + ). The existence of the 
H was predicted by Jaffe in 1977 [1] in the framework 
of the quark-bag model. Its mass was originally esti- 
mated to be around 2150 MeV, making it stable toward 
strong decay to two A particles. Since then, there have 
been many theoretical efforts to determine its mass and 
production cross section and, on the experimental side, 
many inconclusive or unsuccessful attempts to produce 
and detect it. 

Our work is prompted by the possibility that the H is 
lighter than two nucleons. This is motivated by several 
lines of reasoning based on hadron phenomenology and 
non-perturbative QCD modcling[2]. The interpretation 
of A(1405) and A(1520) as bound states of gluon plus uds 
quarks in a flavor-singlet color-octet state, suggests the H 
has properties similar to a glueball and tub rj 1.3 — 1.8 
GeV[2, 3]. Adding current quark masses to the origi- 
nal Skyrme model calculation with massless quarks [2, 4] 
gives ttih as 1.8 GcV. Finally, an instanton- liquid calcu- 
lation gives m H = 1780 MeV [5]. 

Being tightly bound, the H is expected to be a spatially 
compact state. Analogy with the glueball and instanton- 
liquid results suggest rn ~ tq w (1/3 — 1/2) r n w 
(1/6 — 1/4) rjy [2, 6]. In the absence of an unquenched, 
high-resolution lattice QCD calculation capable of a re- 
liable determination of the H mass and size, we will take 
t hI t n = 1// with / treated as a parameter expected to 
be in the range 4-6. For a more detailed discussion of 
motivation and properties of a stable H, and a review of 
experimental constraints on such an H, see ref. [2]. 

In this paper we focus on two types of experimental 
constraints on the transition of two baryons to an H in 
a nucleus, Abb - ► A' H X. Experiments observing single 
A decays from double A hypernuclei A A \{7, 8] indicate 
that t(A aa -> A' h X)<t a = 3 lO" 10 sec. In addition, 
if the H is lighter than two nucleons, nuclei are unstable 
toward AS = —2 weak decays producing the H particle 
in the final state. To estimate the rates for these pro- 



cesses requires calculating the overlap of initial and final 
quark wavefunctions in a nucleus. We will model that 
overlap using an Isgur-Karl harmonic oscillator model for 
the baryons and H, and the Bethe-Goldstone wavefunc- 
tion for a nucleus. Our results will be expressed in terms 
of the parameter / = r^/rn and the nuclear hard core 
radius. We will find that the stability of nuclei is the 
more stringent constraint on the properties of a stable 
H, but is acceptably small if the H is sufficiently com- 
pact: rn < 1/4 rjv depending on mass and nuclear hard 
core radius. Adequate suppression of T(A AA — > A' H X) 
requires ^ 1/3 rjv, whether H is stable or not. Thus 
a byproduct of this investigation is that the conventional 
H with mass Irrifq < tub < 2m A may still be viable in 
spite of the observation of double-A hypernuclei, consis- 
tent with the conclusion of ref. [9]. 

This paper is organized as follows. In section II we 
describe in greater detail the two types of experimental 
constraints on the conversion of baryons to an H in a 
nucleus. In section III we calculate lifetimes for these 
transitions, using a non-relativistic harmonic oscillator 
quark model and a phenomenological treatment of the 
weak interactions. The results are summarized in section 
IV. 



II. EXPERIMENTAL CONSTRAINTS 

A. Stability of nuclei 

There are a number or possible reactions by which 
two nucleons can convert to an H in a nucleus. The 
initial state is most likely to be pn or nn in a relative 
s-wave, because in other cases the Coulomb barrier or 
relative orbital angular momentum suppresses the over- 
lap of the nucleons at short distances which is necessary 
to produce the H. If rtiB % 2rriN — nm^ 1 , the final state 



1 Throughout, we use this shorthand for the more precise inequal- 
ity mu < ttia — ra A i — rax where mx is the minimum invariant 
mass of the final decay products. 
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can be Htt + or Htt° and n — 1 pions with total charge 
0. For m H > 1740 McV, the most important reactions 
are pn — > He + v e or the radiative-doubly-weak reaction 
nn — > i?7. 

The most sensitive experiments to place a limit on the 
stability of nuclei are proton decay experiments. Super 
Kamiokande (SuperK), places the most stringent con- 
straint due to its large mass; it is a water Cerenkov de- 
tector with a 22.5 kiloton fiducial mass, corresponding to 
8 10 32 oxygen nuclei. SuperK is sensitive to proton de- 
cay events in over 40 specific proton decay channels[10]. 
Since the signatures for the transition of two nucleons 
to the H are substantially different from the monitored 
transitions, a specific analysis by SuperK is needed to 
place a limit. We will discuss the order-of-magnitude of 
the limits which can be anticipated. 

Detection is easiest if the H is light enough to be pro- 
duced with a 7r + or 7r°. The efficiency of SuperK to 
detect neutral pions, in the energy range of interest (KE 
= 0-<~ 300 MeV), is around 70 percent. In the case that 
a 7T+ is emitted, it can charge exchange to ir° within 
the detector, or be directly detected as a non-showering 
muon-like particle with similar efficiency. More difficult 
is the most interesting mass range mj?> 1740 MeV, for 
which the dominant channel pn — > He + v gives an elec- 
tron with E ~ (2m N - m H )/2<70 MeV. With a rate 
of order a smaller, the nn — ► Hj channel would give a 
monochromatic photon with energy (2mjv — tur) < 100 
McV. 

We can estimate SuperK 's probable sensitivity as 
follows. The ultimate background comes primar- 
ily from atmospheric neutrino interactions: vN — > 
N'(e,fi), vN -» N'(e,n) + nn, vN -> uN' + nn, 
which has a rate of about 100 kton - yr -1 . Without a 
strikingly distinct signature, it would be difficult to de- 
tect a signal rate significantly smaller than this, which 
would imply SuperK might be able to achieve a sensitiv- 
ity of order ta nn -*a' h x > fewlO 29 yr. Since the H pro- 
duction signature is not more favorable than the signa- 
tures for proton decay, the SuperK limit on ta nn -*a' h x 
can at best be 0.1r p , where 0.1 is the ratio of Oxygen 
nuclei to protons in water. Detailed study of the spec- 
trum of the background is needed to make a more precise 
statement. We can get a lower limit on the SuperK life- 
time limit by noting that the SuperK trigger rate is a few 
Hz[10], putting an immediate limit to^h+x ^ fewlO 25 
yr, assuming the decays trigger SuperK. 

While SuperK limits depend on specific decay chan- 
nels, three other experiments potentially establish limits 
on the proton lifetime which are independent of the de- 
cay channel[ll]. They place weaker constraints on the 
lifetime, due to their smaller size, but are of interest be- 
cause they measure the stability of nuclei directly. The 
experiments of Dix et al.[12] and Evans et al.[13] are not 
sensitive to two nucleon transitions and thus are not ap- 
plicable to nuclei disintegrating with the emission of an 
H. 

The experiment of Flerov et. al.[14] could in princi- 



ple be sensitive to such transitions. It searched for decay 
products from Th 232 , above the Th natural decay mode 
background of 4.7 MeV a particles, emitted with the rate 
T Q = 0.7 10 _10 yr _1 . The conversion of two nucleons in 
Th 232 could result in the following decay chains: 

Th 232 -» R a f°X; Ra 230 ^!^ Ac 230 + /?(0.99MeV) 
Th 232 - Aci°X; Ac 230 Th 230 + /3(2.7MeV) 

Th 23 „ 2 -> Th 2 ff 30 V; Th2 30 7 ^T Ra226 + a (A.7MeV). 
However in general the transitional nuclear state, de- 
noted e.g., Ra 2 ! , would have additional more compli- 
cated decay chains through excited states. Note that the 
H does not bind to nuclei[15]; it simply recoils with some 
momentum imparted in its production. The Flerov et 
al[14] experiment must have cuts to remove the severe 
background of 4.7 MeV a's. If these cuts do not remove 
the events with production of an H, it would imply the 
limit T Th 232^ H+x > 10 21 yr. Unfortunately ref. [14] does 
not discuss these cuts or the experimental sensitivity in 
detail. An attempt to correspond with the experimental 
group, to determine whether their results are applicable 
to the H, was unsuccessful. 



B. Double A hyper- nuclei detection 

There are five experiments which have reported posi- 
tive results in the search for single A decays from dou- 
ble A hypernuclei. We will describe them briefly The 
three early emulsion based experiments [16-18] suffer 
from ambiguities in the particle identification, and there- 
fore are considered less reliable. In the latest emulsion 
experiment at KEK [8] , a double hypernucleus event has 
been observed and interpreted with good confidence as 
the sequential decay of He AA emitted from a S~ hy- 
peron nuclear capture at rest. The binding energy of 
the double A system is obtained in this experiment to 
be Baa = 1-01 ± 0.2 MeV, in significant disagreement 
with the results of previous emulsion experiments, find- 
ing Baa ~ 4.5 MeV. 

The synchrotron based experiment [7] used the 
(K~ , K + ) reaction on a Be 9 target to produce S=-2 nu- 
clei. That experiment detected pion pairs, coming from 
the same vertex in the Be target. Each pion in a pair in- 
dicates one unit of strangeness change from the (presum- 
ably) di-A system. Peaks in the two pion spectrum have 
been observed, interpreted as corresponding to two kinds 
of decay events. The pion kinetic energies in those peaks 
are (114,133) MeV and (104,114) McV. The first peak 
can be understood as two independent single A decays 
from A A nuclei. The energies of the second peak do not 
correspond to known single A decay energies in hyper- 
nuclei of interest. The proposed explanation [7] is that 
they are pions from the decay of the double A system, 
through some specific He resonance. The required reso- 
nance has not yet been observed experimentally, but its 
existence is considered plausible. This experiment does 
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not suffer from low statistics or inherent ambiguities, and 
one of the measured peaks in the two pion spectrum sug- 
gests observation of consecutive weak decays of a double 
A hyper- nucleus. The binding energy of the double A 
system Baa could not be determined in this experiment. 

The KEK and BNL experiments demonstrate quite 
conclusively in two different techniques, the observa- 
tion of A decays from double A hypernuclei. Therefore 
T Aw^A' H x cannot be much less than w 10~ 10 s. (To give 
a more precise limit on ta aa ^a' h x requires a detailed 
analysis by the experimental teams, taking into account 
the number of hypernuclei produced, the number of ob- 
served A decays, the acceptance, and so on.) As will be 
seen below, this constraint is readily satisfied if the H is 
compact: r\ff<l/3 rjv- 



III. BB TO H TRANSITION RATES IN NUCLEI 

As discussed in the introduction, the H may be consid- 
erably more compact than a nucleon which would sup- 
press its production from a two baryon initial state, due 
to the small overlap of the initial and final wave functions 
in position space. We will estimate this suppression us- 
ing the non-relativistic harmonic oscillator quark model. 
Additional suppression comes from the hard core repul- 
sion in the nucleon- nucleon wave function in the nucleus. 
To take that into account, we use the Bcthe-Goldstone 
relative wave function of two baryons in a nucleus. 

The matrix element for the transition Ann — ► A' H X 
is calculated in the AA pole approximation, as a prod- 
uct of matrix elements for two subprocesses: a transition 
matrix element for formation of the H from a AA system, 
|-M|aa->./? x, times the amplitude for a weak doubly- 
strangeness-changing transition, |.<M|Ariv->AA- The sup- 
pression in the spatial wavefunction overlap enters the 
AA — > H transition. We calculate this part of the transi- 
tion amplitude first. The estimate of \M\nn^aax based 
on weak interaction phenomenology is given afterwords. 



A. Calculation of JA'IIaa^h 



where we have neglected constituent quark mass differ- 
ences. The wave function of baryons can then be written 
in terms of the relative positions of quarks, while the 
center of mass motion is factored out. The relative wave 
function in this model is[19, 20] 



r 2 , f 3 ) = N B exp 



a 



(2) 



where Nb is the normalization factor, o-b = 



V< r l> 



V3Km, and < r 2 B > is the baryon mean charge radius 
squared. Changing variables to 



P = — , A = 



V2 



V6 



(3) 



reduces the wave function to two independent harmonic 
oscillators. In the ground state 



*s(p,A) 



OLB_ V 



exp 



-^(P 2 + A 2 ) 



(4) 



One of the deficiencies of the IK model is that the 
value of the as parameter needed to reproduce the 
mass splittings of lowest lying i + and | + baryons corre- 
sponds to a mean charge radius squared for the proton of 
< r 2 ch >= -±r = 0.49 fm. This is distinctly smaller than 

the experimental value of 0.86 fm. Our results depend 
strongly on the choice of as and therefore we should keep 
in mind this problem. Another concern is the applicabil- 
ity of the non-relativistic IK model in describing quark 
systems, especially in the case of the tightly bound H. 
With th/tn — 1//, the quark momenta in the H are 
w / times higher than in the nucleon, which makes the 
non-relativistic approach more questionable than in the 
case of nucleons. Nevertheless we adopt the IK model 
because it offers a tractable way of obtaining a qualita- 
tive estimate of the effect of the small size of the H on the 
transition rate, and there is no other alternative available 
at this time. 

We fix the wave function for the H particle starting 
from the same Hamiltonian (1), but generalized to a six 
quark system. For the relative motion part this gives 



We calculate | M. \ aa-^h in position space as the overlap 
of the H and AA wave functions in the Isgur-Karl (IK) 
non-relativistic harmonic oscillator quark model[19, 20]. 
We take the A spatial wavefunction to be the same as the 
nuclcon's. For now we are concerned with the dynamics 
of the process and we defer discussion of the suppression 
from the spin-flavor part of the transition amplitude. 

The IK model was designed to reproduce the masses of 
the observed resonances and it has proved to be successful 
in calculating baryon decay rates [19]. In the IK model, 
the quarks in a baryon are described by the Hamiltonian 



-pl) + -K^ <1 (r l -r J f (1) 



= Nh exp 



,2 6 



i<j 



(5) 



The space part of the matrix element < A' h \Aaa > is 
given by the integral 



17 rf 3 n^(l, 2, 3)* A (4, 5, 6)* ff (l, 2, 3, 4, 5, 6). (6) 

i=i 



We can rewrite this in a more convenient form, changing 
variables to 

n,f2,f3,f 4 ,f 5 ,f 6 -> p a ,X a 7 p*,X b ,a,R C M (7) 
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where p* 2 ^ and A"^- 1 are defined as in eq (3), with a(b) 
referring to coordinates 1, 2, 3 (4, 5, 6). (Since we are ig- 
noring the flavor-spin part of the wavefunction, we can 
consider the six quarks as distinguishable and not worry 
about fermi statistics at this stage.) We also define the 
center-of-mass position and the separation, a, between 
initial baryons a and 6: 



Rc 



A I 



n CM + U CM -> _ pa p6 

- , a — n CM — n CM . 



(8) 



Using these variables, the H ground state wave function 
becomes 



/ o \ 1/4 / \ 15/2 

•--(*) (r. 

, a 2 H .^2 -2 -2 -.2 3 

x exp[--^-(p a +A a + p b +X b +2") 



(9) 



and the overlap of the space wave functions is given by 

\M\AA^H= [ J] ^V^ 3 AM 3 a i>A VW (10) 



-a,b 



where the center of mass dependence has been factored 

out, Va' 6 = 4>A b (p a ' b A a ' b ), and Vw = Vw(a) is the 
relative wavefunction function of the two A's in the nu- 
cleus. The integration over the center of mass position of 
the system gives a 3 dimensional momentum delta func- 
tion. In the case of pion or lcpton emission, plane waves 
of the emitted particles should be included in the inte- 
grand. For brevity we use here the zero momentum trans- 
fer, k = approximation, which we have checked holds 
with good accuracy; this is not surprising since typical 
momenta are <0.3 GeV. 

To describe two A's or nucleons in a nucleus we will 
use solutions of the Bruecker-Bethe-Goldston equation 
describing the interaction of a pair of fermions in an in- 
dependent pair approximation; see, e.g., [21]. The two 
particle potential in a nucleus is poorly known at short 
distances. Measurements (the observed deuteron form 
factors, the sums of longitudinal response of light nu- 
clei,...) only constrain the two-nucleon potentials and 
the wave functions they predict at internucleon distances 
larger than 0.7 fm [22]. The Bethe-Goldstone equation 
can be solved analytically when a hard-core potential is 
used. While the hard-core form is surely only approxi- 
mate, it is useful for our purposes because it enables us to 
isolate the sensitivity of the results to the short-distance 
behavior of the wavefunction. We stress again, that more 
"realistic" wavefunctions arc in fact not experimentally 
constrained for distances below 0.7 fm. Rather, their 
form at short distance is chosen for technical convenience 
or aesthetics. 

For the s-wave, the B-G wavefunction is 



*bg(«) = 




a > 
a < 



(11) 



where is the hard core radius. The function u vanishes 
at the hard core surface by construction. It then rapidly 
approaches the unperturbed value 1, crossing over that 
value at the so called "healing distance" . Expressions for 
u and Nbg can be found in [21]. 

After performing the Gaussian integrals analytically, 
the overlap of the space wave functions becomes 



\M\aa^h 



V 



1 + /2 

/OG 



1/4 



v^5F 



3/2 



(12) 



3 ^u(k F a) _s n 2 



kpa 



where the factor 1 /4 comes from the probability that two 
nucleons are in a relative s-wave, and / is the previously- 
introduced ratio of nucleon to H radius; olb = / ag. 
Since Nbg has dimensions V^ 1 ^ 2 the spatial overlap 
I-MIaa-*// is a dimcnsionless quantity, characterized by 
the ratio /, the Isgur-Karl oscillator parameter as, and 
the value of the hard core radius. It is shown in Fig. 1 
for a range values of hard-core radius and /, using the 
standard value of as for the IK model[20]. 
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Figure 1: Logio of I^Iaa->.w versus hard core radius in fm, 
for ratio / = Rm/Rh = 4, 5, 6, 7. 



B. Lifetime of doubly-strange nuclei 

We can now estimate the decay rate of a doubly- 
strange nucleus: 



K 2 (2ttY 



in,, 



2(2m AA ) 



(13) 



X ®2\M\Ia^h- 



where $2 is the two body phase final space factor, defined 
as in [11], and toaa is the invariant mass of the A's, 
s=s 2mA- The factor K contains the transition element in 
spin flavor space. It can be estimated by counting the 
total number of flavor-spin states a uuddss system can 
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occupy, and taking K 2 to be the fraction of those states 
which has the correct quantum numbers to form the H. 
That gives K 2 — 1/1440. Thus we obtain the lifetime 
estimate 



1.4 



K 2 \M\ 2 



-10 



-21 



(14) 
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Figure 3: Some relevant weak transitions for NN — > HX 



where the phase space factor was calculated for inn = 1.5 
GcV. 

Fig. 2 shows \M\\ A _>jf in the range of / and hard- 
core radius where its value is in the neighborhood of the 



experimental limits. Evidently, \M\ 



AA~>H 



< 10~ 8 is sat- 
isfied even for relatively large H, e.g., r#<l/3 for 
the canonical choice 0.4 fm for hard-core radius. This 
suppresses T(Aa\ — > A' H X) sufficiently that some A's in 
a doublc-A hypernucleus will decay prior to formation of 
an H. Thus the observation of single A decay products 
from double-A hypernuclei cannot be taken to exclude 
the existence of an H with mass below 2mA unless it can 
be demonstrated that r# > 1/3 r^. 
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Figure 2: Logio of \M\\ A ^ H versus hard core radius in fm, 
for 1=2, 3, 4. 



C. Calculation of the \M\bb^aax matrix element 

Transition of a two nucleon system to AA requires two 
strangeness changing weak reactions. Possible AS = 1 
sub-processes to consider are a weak transition with 
emission of a pion or lepton pair and an internal weak 
transition. These are illustrated in Fig. 3 for a three 
quark system. We estimate the amplitude for each of the 
sub-processes and calculate the overall matrix element 
for transition to the AA system as a product of the 
sub-process amplitudes. 



The matrix element for weak pion emission is esti- 
mated from the A — > Ntt rate: 



\M\ 



1 2m A 1 



A^Nt. 



(2tt) 4 $ 2 t\ — >Nn 



= 0.8 x 10~ 12 GeV 2 . 

(15) 

By crossing symmetry this is equal to the de- 
sired \-M-\%^a-r' in the approximation of momentum- 
independence which should be valid for the small mo- 
menta in this application. Analogously, for lepton pair 
emission we have 



,..|2 1 2m A 1 

\ JVl \\->Nev — TTTTi If, Z 

(2irp <P 3 TA — >Nev 



= 3.0 X 10~ 12 . (16) 



The matrix element for internal conversion, (uds) — > 
(udd), is proportional to the spatial nucleon wave func- 
tion when two quarks are at the same point: 



|A4|a^jv ~< M s (ri-r 2 )\ip N > 



Gf sin C cos 6 C 



m„ 



(17) 



where m q is the quark mass introduced in order to make 
the 4 point vertex amplitude dimcnsionlcss[23]. The ex- 
pectation value of the delta function can be calculated in 
the harmonic oscillator model to be 



< ipA\S 3 {fi-f 2 )\ip N > 



1 = 0.4xl0- 2 GeV 3 . 

(18) 

The delta function term can be also inferred phenomeno- 
logically in the following way, as suggested in [23]. The 
Fermi spin-spin interaction has a contact character de- 
pending on <fi<J2/rn 2 S(fi — F2), and therefore the delta 
function matrix element can be determined in terms of 
electromagnetic or strong hyperfine splitting: 

(m s o - m s +) - (m„ - m p ) = a— < S 3 (f 1 - r 2 ) >(19) 

m A -m N = a S^2 < ^(^i ~ > ( 20 ) 

where m q = taken to be toat/3 is the quark mass. Using 
the first form to avoid the issue of scale dependence of 
as leads to a value three times larger than predicted by 
the method used in eqn 18, namely: 

< ipA\S 3 (n - r 2 )\^N > = 1.2 x 10~ 2 GeV 3 . (21) 

We average of the expectation values (18) and (21) and 
find 



\M\ 2 A ^ N 



4.4 x 10~ 15 . 



(22) 
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In this way we have roughly estimated all the matrix 
elements for the relevant sub-processes based on weak- 
interaction phenomenology. 



D. Nuclear decay rates 

NN — ► HX requires two weak reactions. For the process 
Ann — ► A' h ttit, the rate is thus approximately 



-T A N N ^A' H 1T7T 



K 



2 (2^) 4 



(23) 



2(2m N ) 

x f \M\l^\M\ A ^ H 
\ (2mA - mn) 2 

where the denominator is introduced in the spirit of the 
AA pole approximation, to make the 4 point vertex am- 
plitude dimensionless. The lifetime for this decay is 



0.03 



K*\M\ 2 



(24) 



AA^H 



taking mn = 1.5 GeV in the phase space factor. For the 
process with one pion emission and an internal conver- 
sion, the rate estimate is 



2(2m N ) 

x {\M\ n ^\M\ n ^a\M\aa^h) 2 
leading to the lifetime for m# = 1.5 GeV of 

2 x 10~ 3 



(25) 



+A'tt 



K 2 \M\l A ^ H 



yr. 



(26) 



If niH ^ 1740 MeV, pion emission in a nucleus is kinc- 
matically suppressed and the relevant final states are e + v 
or 7; we now calculate these rates, taking m H = 1.8 GeV. 
For the transition Ann — * A' H ev, the rate is 



Fa* 



*A'„ev 



K 2 



(27T) 4 



.a. 



2{2m N ) 

X (\M\N^Ae^\M\N^A\M\ A A^ H ) 2 ■ 

In this case, the nuclear lifetime is 

70 



(27) 



T~Ann— *A'„ev 



yr. 



K 2 \M\\ A ^ H 
For Ann — ► ^h1^ the rate is approximated as 



^A NN ^A' H 7 



K 2 (2ir 



A a E Mm q 
2{2m N ) 



(28) 



(29) 



x <S> 2 (\M\ 2 n^a\M\aa^h) 2 
leading to the lifetime 

2 10 3 

ta nn ^a> h1 « K *\M\l A ^ H 



(30) 



One sees from Fig. 1 that a lifetime bound of > few 10 29 
yr is not a very stringent constraint on this scenario if 
mn is large enough that pion final states are not al- 
lowed. E.g., with K 2 = 1/1440 the rhs of cqn (28) is 
> few 10 29 yr, for the a hard core radius of 0.45 fm and 
Th ~ 1/5 vn - in the middle of the range expected based 
on the glueball analogy. If thh is light enough to per- 
mit pion production, experimental constraints are much 
more powerful, mn < 1740 MeV is disfavored but not 
excluded; the allowed region in the /-hard core radius 
plane may be reasonable, depending on how strong lim- 
its SuperK can give. 



TABLE I: The final particles and momenta for nucleon- 
nucleon transitions to H in nuclei. For the 3-body final states 
marked with *, the momentum given is for the configuration 
with H produced at rest. 



mass 


final state 


final momenta 


partial lifetime 


m H [GeV] 


A' + H + 


p [MeV] 


xK 2 \M\l K ^ H [yr] 


1.5 


7T 


318 


2 1CT 3 


1.5 


7T7T 


170* 


0.03 


1.8 


eu 


48* 


70 


1.8 


7 


96 


2 10 3 



IV. CONCLUSIONS 

We have considered the stability of nuclei and hyper- 
nuclei with respect to conversion to an H dibaryon. If 
the binding of the H dibaryon is strong, possibly result- 
ing in ran < 2tun as conjectured in refs. [2, 3], then the 
size of the H is expected to be much smaller than the 
size of a nucleon and comparable to the size of a glue- 
ball: ru ~tq w (1/6— 1/4) tn- We used the Isgur-Karl 
wavefunctions for quarks in baryons and the H, and the 
Bethe-Goldstone wavefunction for nucleons in a nucleus, 
to obtain a rough estimate of the wavefunction overlap 
for the process Abb — ► A' H X. We find that observation 
of A decays in double-A hypernuclei does not exclude an 
H - stable or not - as long as rn < 1/3 tn- 

Combining our wavefunction overlap estimates with 
phenomcnological weak interaction matrix elements, per- 
mits the lifetime for conversion of nuclei to H to be es- 
timated. These estimates have uncertainties of greater 
than an order of magnitude: the weak interaction ma- 
trix elements are uncertain to a factor of a few, factors 
of order 1 were ignored, a crude statistical estimate for 
the flavor-spin overlap was used, mass scales were set 
to totv/3, and most importantly, the calculation of the 
wavefunction overlap used models which surely oversim- 
plify the physics. While the overlap is highly uncer- 
tain because it depends on nuclear wavefunctions and 
hadronic dynamics which are not adequately understood 
at present, the enormous suppression of H production 
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which wc found in this calculation forces us to conclude 
that an absolutely stable H is not excluded by these con- 
siderations. 

SuperK can place important constraints on the conjec- 
ture of an absolutely stable H, or conceivably discover ev- 
idence of its existence, through observation of the pion(s), 
positron, or photon produced when two nucleons in an 
oxygen nucleus convert to an H. We estimated that Su- 
perK could achieve a lifetime limit r > few 10 29 yr. Until 
the properties of the H and the dynamics of production 
of the H in nuclei are better understood, this limit would 



be insufficient to rule out a stable H. However such a 
sensitivity would access the estimated lifetime range for 
ran > 1740 MeV and rn ~ 1/5 rjv, and an experimental 
search is warranted. 
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